Phonon-induced Exciton Dephasing in Quantum Dot Molecules 
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A new microscopic approach to the optical transitions in quantum dots and quantum dot 
molecules, which accounts for both diagonal and non-diagonal exciton-phonon interaction, is de- 
veloped. The cumulant expansion of the linear polarization is generalized to a multilevel system 
and is applied to calculation of the full time dependence of the polarization and the absorption 
spectrum. In particular, the broadening of zero-phonon lines is evaluated directly. It is found that 
in some range of the dot distance real phonon-assisted transitions between exciton states dominate 
the dephasing, while virtual transitions are of minor importance. The influence of Coulomb inter- 
action, tunneling, and structural asymmetry on the exciton dephasing in quantum dot molecules is 
analyzed. 

PACS numbers: 78.67.Hc, 63.22.+m, 42.50.Md 



Semiconductor quantum dots (QDs) have been con- 
sidered for several years as promising candidates to form 
an elementary building block (qubit) for quantum com- 
puting. Recently quantum dot molecules (QDMs), i.e. 
systems of two quantum-mechanically coupled QDs, have 
been proposed for realization of optically driven quantum 
gates involving two-qubit operations 0, 0, H] • It turns 
out, however, that contrary to their atomic counterparts, 
QDs have a strong temperature-dependent dephasing of 
the optical polarization. Such a decoherence caused by 
the interaction of the electrons with the lattice vibrations 
(phonons) is inevitable in solid state structures thus pre- 
senting a fundamental obstacle for their application in 
quantum computing. 

There has been considerable progress in the under- 
standing of exciton dephasing in QDs after a seminal 
publication by Borri et al. on four-wave mixing mea- 
surements in InGaAs QDs |^. Two important and well 
understood features of the dephasing are: (i) the opti- 
cal polarization experiences a quick initial decay within 
the first few picoseconds after pulsed excitation and (ii) 
at later times it shows a much slower exponential decay. 
In photoluminescence and absorption spectra this mani- 
fests itself as (i) a broadband and (ii) a much narrower 
Lorentzian zero-phonon line (ZPL) with a temperature 
dependent linewidth |f| . Such a behavior of the polar- 
ization is partly described within the widely used inde- 
pendent boson model [|| that allows an analytic solution 
for the case of a single exciton state. It describes satis- 
factorily the broadband (or the initial decay of the po- 
larization). However, in this model there is no long-time 
decay of the polarization (no broadening of the ZPL). 

Recently we have presented a first microscopic calcu- 
lation of the ZPL width in single QDs 0] , taking into ac- 
count virtual phonon-assisted transitions into higher ex- 
citon states and mapping the off-diagonal linear exciton- 
phonon coupling to a diagonal but quadratic Hamilto- 
nian. This is the major mechanism of phonon-induced 



dephasing in single QDs as long as the exciton level dis- 
tance is much larger than the typical energy of the acous- 
tic phonons coupled to the QDs (< 3meV in InGaAs 
QDs). On the contrary, in QDMs the distance between 
the nearest exciton levels can be made arbitrary small if 
the tunnelling between dots is weak enough so that 
the interaction with acoustic phonons can lead to real 
transitions (changing the level occupation). The exper- 
imentally measured excitonic polarization shows a dif- 
ferent behavior in QDMs ||, too: the long-time decay 
is multi-exponential, in contrast to a single-exponential 
one in uncoupled QDs 0. 

To take into account both real and virtual transitions 
we develop in this Letter a new approach for a multilevel 
excitonic system which is coupled to acoustic phonons 
both diagonally and non-diagonally. This allows us to 
calculate the dephasing in QDMs, as well as the full time- 
dependent linear polarization and absorption. Instead of 
the self-energy approach |loj which is more standard in 
;ron-phonon 

0Q. it 

studying the dephasing and, being applied to a multi- 
level system, has to be generalized to a matrix form. 

To calculate the linear polarization we reduce the full 
excitonic basis to the Hilbert space of single exciton 
states |n) (with bare transition energies E n and wave 
functions ^ n )- Then the exciton-phonon Hamiltonian 
takes the form 

(jjqO^a^ + V nm \n)(m\ 

n q nm 



the electron- phonon problem, we use the cumulant ex- 
pansion 0, ll ll| . It is much more advantageous when 
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where a q is the acoustic phonon creation operator, D c ^ 
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is the deformation potential constant of the conduction 
(valence) band, pu is the mass density, u s the sound 
velocity, V the phonon normalization volume, and H = 1. 
The linear polarization is given by 

P(t) = (d(t)d(0)) = <rf ro e- <B "*P„ m (i) (3) 

nm 

where d n = p cv J dr^ n (r,r) are the projections of the 
excitonic dipole moment operator d (|<in| 2 are the oscil- 
lator strengths) . The components of the polarization are 
written as standard perturbation series, 

Pnm{t) = dh dt 2 ... dt 2k 

k=0 Jo Jo Jo 

x e «o„ p ti e ifi pr t 2 e in sm t 2k 

pr ... s 

x (Vn P 0i)Vp r (i 2 ) • ■ ■ V sm (t 2k )), (4) 

where the finite-temperature expectation value is taken 
over the phonon system, and the difference energies are 
^nm E n E rn . The expansion of cxp(-iE n t)P nm (t) 
which is in fact the full exciton Green's function is shown 
diagrammatically in Fig. 1 up to second order, where the 
phonon Green's function (TV n k(t)V pm (t')) (dashed lines) 
depends on four exciton indices. 

Instead of the plain or self-energy summation of dia- 
grams we introduce for each time t the cumulant matrix 
K(t) defined as 

P(t) = e*<*\ (5) 

where P(i) = 1 + pW + P^ + ... is the expansion 
of the polarization matrix P nrn given by Eq.Q). Then 
the corresponding expansion for the cumulant is easily 
generated: 

K(t) = p« + P^ - 1pm pm + . . . ( 6 ) 

Numerically, we restrict ourselves to a finite number of 
exciton levels and diagonalize the cumulant matrix K at 
a given time in order to find the polarization via Eq.JjjJ. 

If all off-diagonal elements of the exciton-phonon in- 
teraction are neglected (M" m = 8 nm M™ n ), the cumulant 
expansion ends already in first order: the contribution of 
all higher terms of the polarization is exactly cancelled 
in the cumulant, Eq.@), by lower order products. This 
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FIG. 1: Diagram representation of the perturbation series 
for the full exciton Green's function up to second order. 



result allows the exact solution of the independent bo- 
son model 0. The inclusion of the non-diagonal inter- 
action leads to non-vanishing terms in the cumulant in 
any order. Still, there is a partial cancellation of dia- 
grams which provides a large time asymptotics of the 
cumulant, K(t) — > — S — iCot — Ft, that is linear in time. 
Consequently, the lineshape of the ZPL is Lorentzian. 
For example, diagrams (1), (2b), and (2c) in Fig. 1 be- 
have linear in time at t 3> L /u s (L is the QD size) , while 
diagram (2a) has a leading t 2 behavior. In the cumulant, 
however, this quadratic term is cancelled exactly by the 
square of diagram (1), P^ P^ jl. 

The broadening of the ZPL (which is absent in the 
independent boson model) is exclusively due to the non- 
diagonal exciton-phonon interaction and appears already 
in first order of the cumulant. Remarkably, the cumu- 
lant expansion reproduces in first order exactly Fermi's 
Golden rule for the real phonon-assisted transitions: 
ri 1} = ^AWc(AP)£ q |MiT<5(A£; - u q ), where the 

ground state dephasing rate r^ 1 is given here for a sys- 
tem with two excitonic levels only. It is simulated for 
a single spherical QD as a function of the level dis- 
tance AE = E 2 - Ex (Fig. 2, dashed curve). As M* 2 
decays with q due to the localization of the exciton 
wave functions (Gauss type in the present model calcula- 
tion), also decays quickly with AE (see the parabola 
in the logarithmic scale). It exhibits a maximum at 
AE = loq ~ u s / L which is a typical energy of phonons 
coupled to the QD. 

In spite of this 'phonon bottleneck' effect, the virtual 
transitions are always present in QDs due to second- 
order diagrams (2b) and (2c) and lead to a non- vanishing 
broadening of the ZPL everywhere. They have been 
taken into account already within our quadratic coupling 
model -7] which is valid in the opposite limit AE ^> ujq 



1 





real transitions 


7 \ " 


— virtual 


present calc. 


A. \ 


















La 





Level distance [meV] 



FIG. 2: Broadening of the ground state ZPL as a function 
of the exciton level distance Ei — E\ in InAs spherical QD, 
calculated with account for real transitions (in first order), 
virtual transitions (according to Ref. 0]) and for both real 
and virtual transitions up to second order in the cumulant. 
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(dash-dotted curve in Fig. 2). In the present calculation 
we do much better than in Ref. Q : We account for both 
real and virtual transitions on an equal footing (up to 
second order in the cumulant) and cover the full range 
of possible values of exciton level distances (Fig. 2, full 
curve). 

To describe excitonic states in a QDM, we restrict our- 
selves in the present calculation to a four-level model. 
We take into account for electron and hole the two low- 
est localized states each. Without Coulomb interaction 
the electron-hole pair state is a direct product of the 
one-particle states which form our basis of four states. 
Then we include the Coulomb interaction and diagonal- 
ize a four-by-four Hamiltonian. Such a four-level model 
is valid as far as the Coulomb matrix elements are smaller 
than the energetic distances to higher confined (p-shell) 
or wetting layer states. 

Since a strictly symmetric QDM would lead by degen- 
eracy to special features 0, we concentrate here on a 
slightly asymmetric situation: the confining potentials of 
the left dot are 2% deeper than those of the right one. 
This is close to the realistic situation, when the In con- 
centration fluctuates from dot to dot. At the same time, 
the shape fluctuations of the QDM are less important for 
shallow QDs studied in the experiments Thus we 

assume both dots to have the same cylindrical form with 
height L z — 1 nm (in the growth direction) adjusted from 
the comparison with experimentally measured transition 
energies 0,0 (taking 92% of In concentration). Given 
that s-,p-, d-, and /-shells in the luminescence spectra of 
QDMs have nearly equidistant positions ^3] , the in-plain 
confining potentials are taken parabolic with Gaussian lo- 
calization lengths of carriers adjusted to l e = 6.0 nm and 
lh = 6.5 nm |l3|. The electronic band parameters are 
taken from Ref. 14| and the acoustic phonon parameters 
are the same as used previously Q- 

While in single dots the Coulomb interaction results 
in a small correction to the polarization decay, in QDMs 
the exciton wave functions (and consequently M™ m ) are 
strongly affected by the Coulomb energy Even more 
important is the influence of the Coulomb interaction 
and asymmetry of the QDM on the exciton transitions 
energies shown in Fig. 3(a) in dependence on (center-to- 
center) dot distance d. 

At short distances d the tunnelling exceeds both the 
Coulomb energy and the asymmetry, and the exciton 
states are well described in terms of one-particle states. 
Like in symmetric QDMs 0, optically active states 
1) = | SS) and |4) = \AA) are formed from, respec- 
tively, symmetric and antisymmetric electron and hole 
states, while the other two, |2) = \SA) and |3) = \AS), 
remain dark. As d increases, the QDM asymmetry and 
the Coulomb interaction mix these symmetric combina- 
tions as is clearly seen from Fig. 3(b). Finally, in the limit 
of large d the two QDs become isolated (no tunnelling) 
and bright states are |1) = \LL) and |2) = \RR) formed 
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FIG. 3: Exciton energies (a) and oscillator strength (b) cal- 
culated within the four-level model of InGaAs QDM with ac- 
count for the Coulomb interaction and 2% of asymmetry, (c) : 
Exciton dephasing rates of the QDM calculated at T — 10 K. 
The meaning of Ec, A e and A^ are given in the text. 



from electrons and holes both localized on the left and on 
the right dot, respectively. The energy splitting between 
them is A e + Ah, where A e (Ah) is the electron (hole) 
asymmetric splitting due to the slight difference between 
the QDs. The other two, |3) and |4), are spatially in- 
direct exciton states which are dark and split off by the 
Coulomb energy Ec- In contrast, in a symmetric QDM 
the two bright states would be separated by the Coulomb 
energy 

The full linear polarization of a QDM is calculated up 
to second order in the cumulant expansion, using the 
described four-level excitonic model. It's Fourier trans- 
form, i.e., the absorption, shown in Fig. 4 contains four 
finite-width Lorentzian lines on the top of broadbands. 
The width of the broadband is of the order of the typ- 
ical energy of phonons participating in the transitions, 
uj ~ 2 meV. Thus, if two levels come close to each other 
and the broadbands start to overlap, the ZPLs get con- 
siderably wider, due to the real phonon-assisted transi- 
tions between neighboring levels. In QDMs this impor- 
tant mechanism of the dephasing is controlled by the tun- 
nelling which induces a level repulsion at short interdot 
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FIG. 4: Absorption spectrum (linear scale) of an asymmetric 
InGaAs QDM calculated at T = 10 K for different interdot 
distances d. The peaks of the ZPLs are truncated. 

distances. 

However, there is another effect which leads to a quite 
unexpected result: Coulomb anticrossing. As d increases, 
the two higher levels come close to each other and should 
exchange phonons more efficiently. Nevertheless, when 
the anticrossing is reached at around d = 8nm (Fig. 4), 
the ZPL width suddenly drops and never restores at 
larger d. This is due to a change of the symmetry of 
states, owing to the Coulomb interaction. At d > 8nm 
states |3) and |4) become more like \LR) and \RL), re- 
spectively [See Fig. 3(b)], and the exciton-phonon matrix 
element between them vanishes by symmetry. Thus, real 
transitions between states |3) and |4) are not allowed any 
more and their dephasing is only due to virtual transi- 
tions into states |1) and |2). 

Dephasing results for all four states are summarized 
in Fig. 3(c). Apart from the features already discussed, 
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FIG. 5: Temperature dependence of the ZPL widths r„ for 
the lowest four exciton levels in d = 6nm QDM. Inset: real 
(virtual) phonon-assisted transitions between exciton states 
shown schematically by full (dashed) arrows. 



there are also oscillations in T clearly seen for levels |2) 
and |4) in the region between d — 4 nm and 8 nm. For an 
explanation, note that the matrix elements M™, Eq.|J2J, 
are Fourier transforms of the electron (hole) probabili- 
ties. When located in different QDs, they carry a factor 
of exp(iqd). Since the typical phonon momentum partic- 
ipating in real transitions is qo ~ l/L z , one could expect 
that r has maxima spaced by a length of order L z . 

The temperature dependence of T n is shown in Fig. 5. 
At d = 6nm, levels |3) and |4) are already close to 
each other, and real phonon-assisted transitions between 
them are possible. As a result, the dephasing rates grow 
quickly with temperature. At the same time, levels |1) 
and 1 2) are far from each other and real transitions are 
suppressed. Still, virtual transitions contribute every- 
where, with no strong dependence on level energies. 

In conclusion, in the present microscopic approach to 
the dephasing in quantum dot molecules, we go beyond 
Fermi's golden rule and quadratic coupling model Q, 
taking into account both real and virtual phonon-assisted 
transitions between exciton levels on equal footing. 
While the dephasing in single quantum dots is mainly 
due to virtual transitions, in quantum dot molecules real 
transitions dominate the dephasing. We show that the 
broadening of the zero-phonon lines calculated for a few 
lowest exciton states depends strongly on interdot dis- 
tance (via tunnelling), electron- hole Coulomb interac- 
tion, and asymmetry of the double-dot potentials. 
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